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Abstract 

We study the Drell-Yan process irN — > fj, + fi~X at large xf using 
perturbative QCD. A higher-twist mechanism suggested by Berger and 
Brodsky is known to qualitatively explain the observed xf dependence 
of the muon angular distribution, but the predicted large xf behavior 
differs quantitatively from observations. We have repeated the model 
calculation taking into account the effects of nonasymptotic kinemat- 
ics. At fixed-target energies we find important corrections which im- 
prove the agreement with data. The asymptotic result of Berger and 
Brodsky is recovered only at much higher energies. We discuss the 
generic reasons for the large corrections at high xp- A proper under- 
standing of the xf — > 1 data would give important information on the 
pion distribution amplitude and exclusive form factor. 
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The Drell-Yan process of inclusive muon pair production in hadronic col- 
lisions |I| is the most important means of determining the valence parton 
content of hadrons other than the nucleon. The kinematics are defined by 
the longitudinal- momentum fraction of the muon pair, Xp, the squared in- 
variant mass of the pair, Q 2 , and the transverse momentum, Q^, of the 
pair. In the Drell-Yan model the reaction is described as the annihilation of 
a massless quark and an antiquark into a virtual photon which decays into 
a muon pair. Because of rotational invariance and parity conservation, the 
polar angle distribution of massless muons in the muon pair rest frame must 
be of the form 

l + Acos 2 #. (1) 

In the Drell-Yan model A = 1 corresponding to the fact that the annihilation 
of on-shell spin-| particles produces a transversely polarized virtual photon. 
This expectation is maintained within small corrections in the presence of 
order a s (Q 2 ) perturbative QCD radiative corrections ||. 

Experimentally, a transverse photon polarization is observed in most of 
the kinematic region. At large xf, however, the photon appears to become 
longitudinally polarized @, £|. In the following, we shall study a higher-twist 
mechanism leading to longitudinal polarization, within a model originally 
suggested by Berger and Brodsky (BB) ||. 

To produce a large xp muon pair, it is necessary that x a , the momentum 
fraction carried by the annihilating projectile parton, is close to 1. Fock 
states F of the projectile H where one parton carries most of the momentum 
have large energy, 

Eh — Ep = AE ~ - kia (2) 

2p(l - x a ) 

where — kj_ a and (l — x a )p are the transverse and longitudinal momenta of the 
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soft parton(s) in the Fock state. In the limit x a — > 1, the lifetime of the Fock 
state, Tp = 1/AE oc 1 — x a is short, and perturbation theory is applicable. In 
the BB model, the leading contribution with a pion projectile is then given 
by the diagrams of Fig. 1, where the spectator valence quark transfers its 
momentum to the active quark via single gluon exchange. Thus the whole 
pion state participates in the scattering, implying longitudinal polarization 
of the virtual photon due to helicity conservation in the pion-virtual photon 
interaction. 

As a first approximation, BB took the nonperturbative wave function 
of the pion to be 5(z — |), i.e., both valence quarks carry half of the pion 
momentum^. In the limit 

z a -l, k 2 ±a /Q 2 ^0 with ki a /Q 2 ~(l-x a ) 2 , (3) 



x a ) 2 (l + cos 2 #) + — ^sin 2 #. (4) 



the differential cross section is proportional to || 

4k_ 2 
~9Q 2 

This gives a simple expression for the parameter A in the angular distribution 

(1 - x a f - 4kj g /9Q 2 

(l-x a ) 2 + 4ki a /9Q 2 ' 1 ' 

It is clear that if (1 — x a ) 2 vanishes faster than k. 2 ±a /Q 2 , the angular distribu- 
tion turns into sin 2 #, and the parameter A approaches —1. Experimentally, 
a drop in A at large xp is observed but the data lie well above the curve 
derived from the asymptotic expression (||). The purpose of this letter is 
to show that there are important nonasymptotic corrections to the limit- 
ing expression (^) at present energies. The BB model evaluated with exact 
kinematics is, in fact, in rather good agreement with the data. 
1 This assumption was relaxed and the BB model generalized in Ref.[[| 



The muon pair production cross section resulting from the model of Fig. 
1, without taking the limit (^) and after integrating over the azimuthal angle 
of the is 

Q 2 da*~ N ^ + ^ x , , ,_ ftN , 2l _, f dx a dx b d 2 k ±b d 2 kj 



dQ 2 d 2 Q±dxpd cos 9 
x [e 2 u fu/N(x b , kxb) + e 2 d f S/N (x b , k ±b 
x5(Q 2 + Qi + 4p 2 



oc (aa s ^(r ± = 0)) 2 |p| J 



1 - x a ) 2 p 2 



St* 



Qn (1 + cos 2 9) + g 00 sin 2 9 



IPI 



kl b + x 2 bV 2 - sJk\ a + {l-x a ) 2 V 2 
xS(x F -x a + x b )5(Q± - k_L a - k_L 6 



(6) 



where p is the pion momentum in the hadron center-of-mass frame, and we 
have expressed the invariants at the 7r-parton level as 



2p ■ p b = 2|p|( v /ki 6 + x 2 p 2 + x b |p|) 
-2p-px = — 2|p 



(7) 

:i-^)Ipd, (8) 



? ±a + (l-x a ) 2 p 2 

neglecting all masses but Q 2 . In terms of s and t, the diagonal density matrix 
elements Qmm = Q-m-m for the production of a virtual photon with spin 
projection M in the Gottfried-Jackson frame are 

en = , , , W JL -^\s 5 + 3sH + 3s A Q 2 + As 3 t 2 + Qs 3 tQ 2 

(s + ty(2Q z — s — ty L 

+ Ash 3 - 2s 2 t 2 Q 2 - 4s 2 Q 6 + 3st A - 10st 3 Q 2 + 8st 2 Q 4 + t 5 

- 5t 4 Q 2 + 8t 3 Q 4 - 4t 2 Q 6 ] , (9) 
8stQ 2 

g °° = -£TW- (10) 

In the following, we will neglect the target parton transverse momentum kj_ b 
in Eq. (|6|). In terms of the density matrix elements above, the polarization 
parameter A of Eq. (p]) then becomes: 

Qn — Qoo 



A 



011 + 000 



(11) 



In Fig. 2 we show A as a function of xf- Our result (solid line) takes 
the exact kinematics into account. We compare this to the BB limit of Eq. 
(D (dashed line) and to the E615 data §]. We fix k 2 ±a = 0.8 GeV 2 , Q = 
4.5 GeV and pbeam = 252 GeV for both calculations, in accordance with the 
data ||. There is a sizeable difference between the solid and dashed curves, 
indicating the importance of nonasymptotic kinematics at present energies. 
Furthermore, the agreement with the data is clearly better using the exact 
kinematics. At this value of k^ a , we find that the general expression (^) 
gives values of A consistent within AA = 0.1 with the asymptotic limit (0) 
for xp < 0.9 when Q 2 and s are scaled by a common factor > 25. Since the 
mean value of k. 2 ±a should increase with Q 2 , the approach to the limit (|5]) will 
in reality be even slower. 

There are generic reasons for the large finite energy corrections at high xp. 
This can be seen, e.g., from the definition of t in (^). For |p| — > oo at fixed x a 
and kj_ a , the leading term of C(p 2 ) in t cancels, giving t ~ — kj_ a /(l — x a ). 
On the other hand, for x a — > 1 at fixed (albeit large) |p| we get the very 
different result t ~ — 2|p||kj_ a |. A similar sensitivity to the order in which 
the limits are taken can be seen in the argument of the first delta function 
in Eq. ([]). This is the constraint which fixes the relation between the BB 
model variables x a , x^ and the measurable quantities xp, Q 2 . In Ref . Q, the 
data was compared to the BB model as a function of x a , using the relation 
x a = \{xp + \J x 2 F + <5 2 /|p| 2 ), which holds only in the asymptotic limit. In 
Fig. 2, we chose the physical quantity xp instead, to avoid ambiguities related 
to model kinematics. 

At present fixed target energies, |p| = 0(10 GeV), hence (1 — x^) |p| = 
(9(1 GeV) ~ (|Q_l|) for xp = 0.9. Thus it is not surprising to find large 
deviations due to finite energy effects. The improved agreement with data 
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obtained in Fig. 2 when such effects are taken into account is of course 
encouraging. However, it is also an indication that one should carefully 
reconsider the applicability of the twist expansion of QCD, on which the BB 
model is based. 

In the usual, leading twist QCD approach, the time scale t i ~ 2p/Q 2 
of the hard interaction (here qq — > fifi) is much shorter than the lifetime of 
the Fock state, t f ~ 2p(l - x a )/k 2 ±a {c.f. Eq. (@)). In the BB model [§], 
the limit (|3]) was taken such that tj/t f = 0(1 — x a ) is vanishing. Hence 
the factorization between the wave function dynamics and that of the hard 
scattering subprocess is still valid. On the other hand, in the limit 

x a ^l, Q 2 ^oo with k 2 ±a /Q 2 ~ (1 - x a ) (12) 

we have tj ~ r F and the dynamics of the subprocess scattering is inseparable 
from that of the Fock state. In this limit the twist expansion breaks down 
and the two new diagrams shown in Fig. 3 contribute at leading order to the 
muon pair production process. 

At finite energies, it is not obvious which (if any) of the asymptotic limits 
(§), ([12]) is more appropriate. Although one may, as we have done here, take 
into account the exact kinematics within a given model, any application of 
perturbative QCD must still depend on an idealized high energy limit. It 
would be worthwhile to systematically study how the effects of finite energy 
corrections may be minimized, since they reflect the inherent uncertainties of 
the approach. In the present case, we did check numerically that including 
the diagrams of Fig. 3 does not significantly change the prediction for A 
given by the solid curve in Fig. 2. 

Muon pair production is a good test case for studying large x F QCD 
dynamics, since the number of diagrams is relatively small, and data is avail- 
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able. This reaction is moreover used to determine the structure function of 
the pion, as well as the shadowing effects of antiquarks in nuclei ||. The 
reliability of these determinations at large xp depend on a proper under- 
standing of the dominant reaction mechanism. In the limits @, (|T2|) the 
production cross section is not given by the projectile structure function but 
rather by the distribution amplitude, i.e., the valence wave function at van- 
ishing transverse distance between the partons ||, as indicated in Eq. (|6|). 
The pion distribution amplitude is in itself of considerable current interest. 
It is not conclusively settled whether the exclusive pion form factor data can, 
at present energies, be described using the asymptotic hard QCD dynamics 
of Ref. |§ or whether the nonperturbative "Feynman" mechanism [|10J is 
more appropriate. A better understanding of large xp muon pair produc- 
tion could help resolve this issue, by yielding information on the distribution 
amplitude. 

While completing this work, we learned of a related study of the muon 



angular distribution at high xf ||1 1|| ■ These authors are considering the effects 
of terms where the transverse momentum of the pair is comparable to 
Q 2 . In the high energy limit where this ratio is nonvanishing, the lifetime 7> 
of the Fock state is actually much shorter than the time scale tj of the hard 
scattering. This illustrates a further possibility of taking limits, not covered 
byEqs. (|) and 

To summarize, we have found that there are general reasons to expect 
large finite energy corrections when applying QCD expressions that are valid 
in xf — > 1 limits to current data. As a case study, we considered the model 
proposed by Berger and Brodsky ||, and found that this model agrees with 
data on the angular distribution of the muon pair only when the kinematics 
are treated exactly. 
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Figure Captions 



Fig. 1. The relevant Feynman graphs for the Berger and Brodsky || 
higher-twist model of Drell-Yan production. 

Fig. 2. The polarization parameter A of Eq. ([!]) as a function of xp. 
Our calculation flTTD with exact kinematics is shown by the solid line, the 
prediction (Q) of the Berger-Brodsky limit (|3|) by the dashed line, and the 
data [§ by the boxes. We used p beam = 252 GeV, Q = 4.5 GeV, k 2 ±a = 0.8 
GeV 2 and k 2 ±b = 0. 

Fig. 3. The two higher-twist diagrams which, together with the diagrams 
of Fig. 1 , contribute at leading order in the limit ( |T^ ) . 
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